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Abstract 

We study, with various methods, the high temperature phase transition for the 
^-component </> 4 theory in the large N limit. Our results fully confirm a previous 
investigation of the problem, for arbitrary N, with the method of the average poten- 
tial. The phase transition is of the second order with an effectively three-dimensional 
critical behaviour. 



1. Introduction 



The iV-component 4 theory is the prototype for investigations concerning the question of 
symmetry restoration at high temperature. Following the original argument of Kirzhnits 
and Linde on the restoration of gauge symmetry J!]], it was considered in all subsequent 
investigations of the problem [^j, £| M. The physical importance of the theory arises from 
the fact that, for iV = 4, it describes the scalar sector of the electroweak standard model 
in the limit of vanishing gauge and Yukawa couplings, and is relevant for the QCD phase 
transition for vanishing up and down quark mass ||. For a negative bare mass term the 
theory can exhibit spontaneous breaking of the O(N) symmetry to 0(N — 1). It was 
demonstrated in that the original symmetry is restored at sufficiently high temper- 

ature. This was achieved through the perturbative calculation of the effective potential 
J7|, and its generalization for non-zero temperature. However, it was realized in the orig- 
inal works that the perturbative expansion breaks down near the critical temperature. 
This is due to the presence of infrared divergencies which destroy the convergence of the 
perturbative series. As a result the details of the phase transition cannot be discussed in 
any finite order in perturbation theory, and even the order of the transition cannot be 
predicted. An amelioration of the situation was achieved in [ffl with the summation of an 
infinite subclass of diagrams (the "daisy" graphs) in the large N limit, but for small N the 
question remained open. Different perturbative calculations led to predictions of a first 
order phase transition |J, while improved treatments reconfirmed the breakdown of the 
perturbative expansion near the critical temperature || . The situation was eventually re- 
solved in |H| through the method of the average potential which employs renormalization 
group ideas JET], [□], |T3| . This method leads to an evolution equation for the dependence of 
the effective potential on a variable infrared cutoff, which provides a good handling of the 
problematic infrared regime by effectively resumming all the relevant classes of diagrams 
of perturbation theory. The high temperature phase transition for the four-dimensional 
(f) 4 theory was shown to be of the second order belonging to the same universality class 
as the phase transition of the three-dimensional 0(iV)-symmetric a- models. A detailed 
quantitative picture was obtained, including the critical exponents for the behaviour in 
the vicinity of the critical temperature. 

The main purpose of this work is to confirm the results of [1C] through the use of 
alternative methods. As it was explained before, such methods should go beyond any finite 
order in perturbation theory. It is for this reason that the large N limit is particularly 
useful. In this limit the functional integral which defines the effective potential can be 
explicitly evaluated. In an alternative approach the Schwinger- Dyson equations of the 
theory can be solved in the critical region. In this way we shall be able to obtain a 
completely satisfactory picture of the phase transition in order to verify the results of 
[ |T0f . We shall compare with the work of ref. [|14]], where the large N limit of the 4 theory 
has also been discussed. Going beyond the leading 1/N result is more difficult and the 
use of the renormalization group becomes indispensable. We shall demonstrate the power 
of the technique employed in JTIJ by using it in order to calculate the critical exponents 
for small N. 
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The paper is organized as follows: In section 2 we first summarize the formalism of the 
average potential for the iV-component 4 theory. Then we use standard large N tech- 
niques in order to obtain a set of equations which permit a straightforward calculation of 
the average potential. As a result we show how the average potential interpolates between 
the classical and the effective potential. Subsequently we obtain an exact renormalization 
group equation for the average potential, which can be used in order to solve for the 
effective potential with the classical potential as a boundary condition. This equation is 
identical to the one employed in [nj in the large N limit. We also generalize our formal- 



ism to the case of non-zero temperature. In section 3 we compute directly the large N 
limit of the effective potential in the critical region (temperatures close to the critical one) 
without using the concept of the average potential. In section 4 an alternative approach 
is followed, which solves the Schwinger-Dyson equations in order to reproduce the results 
of section 3. Sections 3 and 4 are self-contained and can be read independently of the 
rest of the paper. In section 5 we demonstrate the exact agreement of the results of the 
previous two sections with . Then we make use of the renormalization group equation 
in order to go beyond the leading order ml/N and obtain the critical exponents for small 
N. Our conclusions are given in section 6. 

2. The average potential in the large N limit 

We consider a theory of N real scalar fields x a ■, a — 1---N, with an 0(iV)-invariant action 
of the form 

S[ X ] = J d d x{U, X a d»Xa + V(0}, £ = \x a Xa- (2.1) 

For the time being the dimensionality of the Euclidean space-time, d, is kept arbitrary. 
Following ref. [ , we obtain from the above fundamental (or "microscopic") action S an 
average (or "block spin") action 

exp{-r„#]} = J VxP k {<P,x]eM-Six]}- (2.2) 

The constraint operator Pk[<f>, x] enforces the average of x over a volume with size ~ k~ d 
to equal the average field up to small fluctuations. It is given by the expression 

P k [<j>,x) =exp{-i|^x[0 a -/ fe (-n)x a ]P(-n)[0 a -/ fc (-n)Xa]}. (2.3) 

Here fk{l 2 ) is the Fourier transform of the function used in order to implement the 
averaging of x an d different "averaging schemes" correspond to different choices of 
We consider the two parameter family of schemes 

/ fc (g 2 )=exp{-a^ }, (2.4) 



(with a ~ 1 and b > 2 |TT|). P is given by 

P{q 2 )=d 2 [l-fl{q 2 )]-\ (2-5) 
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and has a minimum 

mmP(g 2 ) ~ k 2 . (2.6) 
In eq. ( |2.3|) q 2 is replaced by — □ and the resulting operator acts on x a { x ) i n position 



space. The right-hand-side of eq. ( |2.3| ) contributes a factor q 2 f%(l — f 2 ) 1 to the term 



quadratic in x a - Combining this term with the kinetic term from the classical action ( |2.1 
we obtain 

exp{-I\#]} =exp{-i| d d xrP{-^)(t>a} 

(2.7) 

In order to implement the large N approximation we introduce two (scalar, O(N)- 
singlet) auxiliary fields B(x) and C{x) and write [15 | 

exp{- J d d xV(0} = JvBVCexpl^ - J d d x(C(£- B) + V(B)) j. (2.8) 

Here we have used a rotated integration contour for the delta-functional 

5[X] = J VCexp{ - J d d xCX). (2.9) 
By inserting (|2.8|) into (|2.7[) the action becomes quadratic in x an d the integration can 



be performed easily. We obtain 

exp{-r fe [0]} = J VBVCexp{ - S aux [B, C; 0]}, (2.10) 

with 

S aU x[B, C; = 

| ^{-i<^ ^(-□)p(-D) __J__ /fc (-D)p(-D) _ p(_ D )j 0a + V (B) - SC| 

+ ylndet[P(-D) +C]. (2.11) 
The last term of eq. (|2 . 1 1| ) originates in the internal x-l°°P s an d is proportional to N. 



The other terms in (|2.11| ) are not manifestly proportional to N, but this can be achieved 



by assuming an appropriate N dependence for the fields: ~ VN, B ~ N For a 
polynomial potential of the form 

V(B)=Y,92 P B p , (2.12) 



* Formally one has to rescale B —> B / N and consider S aux [B, C\ v N<f>] instead of S aux [B,C';(j)\. 
Our naive assumption for the N dependence of 4> and B reproduces the effects of this procedure, while 
simplifying the counting of the N dependence of the various terms. 
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the contribution of the potential is ~ iV if the generalized couplings g^ p have an N 
dependence 

92 P -N 1 ~P. (2.13) 

More specifically, the quartic coupling g& = A/2 must scale as g± ~ 1/N for large N [[RJ. 

Since now S aux ~ N, in the limit N —> oc the functional integral ( [2.1 1|) is dominated 
by its stationary point. Therefore 

r k [<i>] = s WIB [B(<i>),c{<f>y,<f>], (2.i4) 

where B and C are the ^-dependent solutions of the saddle point conditions 

6 -^[B,C]=0, (2.15) 

6 4^[B,C}=0. (2.16) 

We are interested in translationally invariant states of the theory which correspond to 
constant 0, B, C . For such states the average action Y k is related to the average potential 
U k through the expression [] 

T h [<f>] = J d d xU k (p). (2.17) 
Here p is an abbreviation which we shall frequently use in the following 

P=\F<l>a. (2-18) 



For constant <fi and C the first term in eq. ( |2.11|) simplifies considerably since 

= pC. (2.19) 



fkPp jjfkP ~ P 



In order to derive the above expression we have exploited that P{q 2 ) diverges P{q 2 ) ~ 
(g 2 ) 1 b for q 2 0. We also notice that for constant C the operator P(—D)+C is diagonal 
in momentum space. Combining all of the above we find 

N r d d a r -i 
U k (B, C; p) = V(B) + (p-B)C+-j A —^- d \n [P(q 2 ) + C] . (2.20) 

An ultraviolet cutoff A is implied for the momentum integration in the previous expression, 
and B, C are assumed to be solutions of the saddle point conditions ( |2.15| ), ( |2.16| ), which 
read 

C=V'(B), (2.21) 

N [ d d q I . 

B =p+— \ , -b , , 2.22 

H 2 J a 2tt ) d P{q 2 ) +V'{B) y 1 



' The wavefunction renormalization is zero in the large N limit. For this reason we do not explicitly 
discuss the derivative terms in the effective action. 
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The last equation gives B as a function of p. Eliminating C from the expression for U k , 
we finally obtain for the average potential in the large N limit 



U k (p) = V(B) + (p - B)V'(B) + *t f \n 



2 Ja (2vr) d 



P{q 2 ) + V\B) 



(2.23) 



with B(p) determined through the "gap equation" ( gjgg ). For k = A eqs. ( gjgg ), (|2^2l|) 
assume the simple form 



f/ A (p) =y(p). 



(2.24) 



In this limit the average potential coincides with the classical potential. For k — > the 
modified inverse propagator P(g 2 ) of eqs. (|2.4|) , ( |2.5|) approaches the standard form 
P(q 2 ) = q 2 . In the same time the average potential U k approaches the convex effective 



potential U [0], |12[. In this limit eqs. ( |2.22j ), ( |2.23| ) coincide with the well known JT5[ 
results for the effective potential in the leading order of the 1/N expansion. An interesting 
formula holds for the p-derivative of Uk- Differentiating ( |2.23j ) and making use of ( |2.22| ) 
we obtain 

j-U k {p) = U' k (p) = V (B(p)) . (2.25) 

We now proceed to derive a renormalization group equation for the average potential 
U k , which is exact in the limit iV —>■ oo. We shall see that this equation coincides with 
the one obtained through the "renormalization group improved" one-loop calculation 



IT] , 12, rjj]. We start by investigating how U k responds to a change of the scale k. Taking 



the logarithmic derivative with respect to k of eq. (|2.23|) and making use of the "gap 
equation" ( |2.22|) we find 



d d q 



A (27r) d P(q 2 ) + V'(B) dk 



k-^Piq 2 )- 



(2.26) 



Combining eqs. ( |2.25| ), ( |2.26| ) leads to the following integro-differential equation 



d d q 1 , d 2 

'a (27r)^P(g 2 ) + U' k (p) dk WJ 



d 

12" 



P(x) + U' k (p) dk 



k4rP(x) 



where 



Vd 




(2.27) 



(2.28) 



This is the desired exact evolution equation for U k . It describes how the average potential 
evolves as the infrared cutoff is lowered from k = A to k = and successive modes are 
integrated out. With the initial condition U k= \(p) = V(p), it provides the means for 
calculating the effective potential U(p) = U k= o(p). The form of the inverse propagator 
P{q 2 ), given by eqs. ( |2.4| ), ( |2.5| ), ( |2.6| ), provides for an effective infrared cutoff for all the 
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X- modes with q 2 < k 2 . Also for modes with q 2 > k 2 the inverse propagator is essentially 
independent of k. As a result, only the modes with q 2 ~ k 2 contribute to the integral 
in ( |2.27[ ). Eq. (|2.27| ) can be compared to the evolution equation obtained through the 
"renormalization group improved" one-loop calculation and used for investigations of the 
average potential in [JXTJ, O, ffO]. The two equations are identical in the limit A" — > oo. 
We finally consider a specific classical potential V(£) 



v(0 



Jl 2 > 0. 



(2.29) 



For the above form of the classical potential (|2.22|) , (|2.23|) , ( |2.25|) result in the following 
equation for the average potential 



U' k (p) = -fi 2 + \p + ^\ 



d d q 



(2.30) 



lA(2<K) d P(q 2 ) + U' k (p) 

For k — 0, where P(q 2 ) = q 2 , this is the well known result for the effective potential in 
leading order in the 1/N expansion |L5| |. Here it is generalized to the average potential 
for arbitrary values of the scale k. By performing a first mass renormalization 

d d q 1 



i4 



we find 



U' k ( P ) 



2 \ -A v 

"Mi + Ap + y A 



d d g q 2 -P(q 2 
A (27r) d ? 2 P(g 2 ) 



a (27r) d g 2 



yA^(p) 



f2.311 



/A(2vr)^)(P(g 2 ) + ^(p))- 

(2.32) 

The integrals in the above expression are ultraviolet and infrared finite in three dimen- 
sions. We can, therefore, take A — ► oo in ( |2.32|) . (In four dimensions additional logarith- 
mic renormalizations must be performed.) Evaluation of the integrals gives an implicit 
equation for U' k (p) as function of p, k, /i 2 and A. This can be used as a check of results 
obtained through different approaches in the limit N —>■ oo. In particular, we find for the 
k dependence of the quartic coupling 




d d q 



1 



-i 



(2.33) 



lA(2ir) d (Ptf) + U k (p)) 2 J 

In order to extend our discussion to the non-zero temperature problem we only need to 
recall that, in Euclidean formalism, non-zero temperature T results in periodic boundary 
conditions in the time direction (for bosonic fields), with periodicity 1/T JE7|. This leads 
to a discrete spectrum for the zero component of the momentum q 



m 



0,±1,±2,... 



(2.34) 



As a consequence the integration over qo is replaced by summation over the discrete 
spectrum 

d d q msr^ f ^ d-1( f 



(27T) 



TJ2 



(2vr) rf - 1 



(2.35) 
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With the above remarks in mind we can easily generalize our results in order to take into 
account the temperature effects. Eq. (|2.30|) is replaced by 

U>{p, T) = -p 2 + Xp + —XT £ j A (27r) ,- lp(g -, + 47r2m2T2) + [/ , (p;T) (2-36) 

in an obvious notation. (The cutoff A is taken A ^> T.) Again, the usual tempera- 
ture dependent effective potential 0-0 obtains from U^{p,T) in the limit k — > 0. The 
generalization of the evolution equation ( |2.27| ) reads 

d TT ^ f d d - l q 1 x^hj , „ 2 2^ 



k dk Uk ^ T) = 2 T ? A (2^ P(g^ + + Ufa, T) k Tk Pif + 47rWT 

(2.37) 

The results of this section can be summarized in eqs. ( |2.22|) , (|2.23| ) and ( 2.27]) . The 



first two determine the average potential Uk(p) for any value of the scale k in the large N 
limit. The third one provides an exact equation for the evolution of the average potential 
as a function of k in the same limit. We also point out the simple relation between the 
derivative of the average potential and the classical one ( |2.25| ). In the specific case of 
the classical potential ( |2.29| ), eqs. ( |2.22| ), ( |2.23| ) are replaced by ( |2.30| ). For non-zero 



temperature eqs. ( |2.30D , Q2.27p take the form ( |2.36D , fl2.37p respectively. In the following 
sections we shall use these results in order to study the high temperature phase transition 
for the four-dimensional, 0(A r )-symmetric 4 theory. 



3. The high temperature phase transition in the large 
N limit 

From this point on we restrict our discussion to four dimensions and to the classical 
potential ( |2 .29| ) . We first consider eq. ( |2 .3U| ) for k — 0. In this limit the average 



potential Uk is equal to the effective potential U and the inverse propagator P(q 2 ) assumes 
its standard form P(q 2 ) = q 2 . We thus obtain the well known result for the effective 



potential in the leading order of the 1/N expansion [15 



^ W = _^ + V + £^_l_, (3.!) 

where we remind the reader that p = |0 a a and the derivatives of U are taken with 
respect to p. The generalization of the above expression for non-zero temperature follows 
2.36| ) and reads 



UUT) = -fi + Ap + f ATg I ^ - _ _ _ . (3.2) 

It is convenient to define the functions 

d 3 q 



/(T) ~ T ? /a (2tt)3 f + ATT 2 m 2 T 2 (3-3) 
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I n {w, T)-TJ2J^ (2 ^ 3 ^ + 4?r2m2T2 + w)n ■ (3.4) 



Then eq. ( |3.2| ) can be written as 
U'( P ,T) 



y A/(T) - /, 2 



+ V+^A[/ 1 (t/'(p ) T),T)-/(T)]. (3.5) 



We are interested in solving the above equation near the phase transition. This correponds 
to a temperature region in which typically U'(p,T) <C T 2 near the minimum of the 
potential. We can, therefore, use a high temperature expansion in terms of ^ for the 
evaluation of I n (w, T) in this regime. In the opposite limit of zero temperature the discrete 
summation can be replaced by an integration over a continuous q° (see ( |2.34j )) and the 



integrals become four-dimensional. Integrals similar to /(T), I n (w, T) have been discussed 
extensively in the literature pi |T7|]. We briefly review here the results which are relevant 



for our investigation. I(T) is given by 

A 2 T 2 

I(T) = T7- 2 + ^- (3-6) 

l07T Z 12 

The first term corresponds to the quadratic divergence of the zero temperature theory, 
while for the evaluation of the second term we have made use of A > T. Notice that 
I(T) receives contributions from all values of m in the infinite sum. In this sense it 
incorporates effects of the effective three-dimensional theory (modes with m = 0) as well 
as the four-dimensional one (m ^ 0). For the difference Ii(w,T) — I(T) we find M 

««, T) - ,m = - JLtVS - ^ in (j^) , (3.7) 

with c = exp(| + ln(47r) — 7) ~ 11.6. For w <C T 2 the above expression gives the first 
two terms of the high temperature expansion. The leading term in the high temperature 
expansion ~ Ty/w comes entirely from the m = term in the infinite sum. As a result eq. 



( |3.7|) in the leading order for small j% reflects only the dynamics of the effective three- 
dimensional theory. In contrast, the logarithmic term is typical for the four-dimensional 
renormalization effects from the momentum range between cT and A. Notice that we 
have included a w-dependence in the argument of the logarithm which is negligible for 
small |j. As a result ( |3.7| ) is also exact for T = 0. (The first term in ( |3.7| ) is not present 
for large but finite but since this term is small in this region it does not affect our 
discussion.) Eq. ( |3.7| ) is therefore valid both for the high and low temperature regime. 
The corrections to ( |3.7[ ) only affect an intermediate "threshold" region where w ~ (cT) 2 . 
The integrals I n (w,T) for n > 1 can be obtained by differentiation of eq. (|3.7|) with 
respect to w. For later use we give the next two integrals 

T , -, s IT 1 , / A 2 \ 1 w 
I 2 {w, T)=—^ + —\n( - 2rrr , 2 , 3.8 

87T y/W 16n z \{cT) z + wJ 16n z (cT) z + w 

y^l^^S (3-9) 



"32tt w I 32vr 2 ((gT) 2 + 



w 
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Again the first term in the above expressions is determined by the effective three-dimensional 
theory. 

Combining eqs. (|3.5[ ), ( |3.6| ) and keeping only the leading contribution in ( |3.7| ) we 
obtain for T 2 > U'{p,T) 



U'(p,T)=NX^ 



Ji 2 -NX 



A 2 



32tt 2 



+ X P -NX^U'(p,T). 



(3.10) 



The ultraviolet divergences can be absorbed in the (T = 0) renormalization of the mass 
parameter and coupling. We define the renormalized parameters | for the zero temperature 
four-dimensional theory |[l(| 



Xr X 
Xr =- 



32tt 2 ' 
A 



NX 
32vr 5 



(3.11) 
(3.12) 



where M is an appropriate renormalization scale M <C A. Eq. ( |3.10| ) now reads 



U'(p,T) = X 



T 2 a 2 
24 X R 



+ Xp-NX^U'(p,T). 



(3.13) 



The critical temperature T cr is defined as the temperature at which U'(0,T cr ) = 0. This 
gives the well known result f|-|| 



T 2 = 24/4 



(3.14) 



We concentrate on temperatures T ~ T cr and /) C T 2 , for which we obtain from eq. 
( |3.13j ) in leading order 



) I 3 T NT, 



(3.15) 



This result justifies a posteriori the high temperature expansion for the integrals I n . 
Notice how any reference to the bare parameters has disappeared in the above expression. 
Simple integration gives for U(p, T) 



U(P,T) = - 



7T 2 T 2 -TV 



T 



P + 



1 8tt 2 T 2 - T c 2 ^ 2 ( 1 64vr 2 1 ^ 3 



A^ 3 T 2 



N 2 3 T 2 



(3.16) 



The last equation gives the effective potential in the critical region T ~ T cr , for large A^ 
and including the p 3 (0 6 ) term. Eq. ( |3.16| ) is in agreement with [|rj| in the leading order 

■"■Note that these parameters correspond to T = and do not coincide with the quadratic and quartic 
term in the temperature dependent effective potential. 
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of the 1/N expansion. It is also completely analogous to results obtained in the context 
of three-dimensional field theory [jTHJ . We observe an N dependence for the renormalized 
couplings in agreement with the general remarks of section 2 (See eq. ( |2.13| ).) Notice that 
the only surviving reference to the parameters of the four-dimensional theory exists in the 
definition of the critical temperature ( |3.14| ). According to our discussion of the integrals 
I, I n , the critical temperature is determined by the dynamics of both the four-dimensional 
and three-dimensional theory. On the contrary, the form of the potential in the critical 
region (and therefore the critical behaviour) reflects only the effective three-dimensional 
theory. 



Eq. (|3.16|) predicts a second order phase transition. For T > T cr the theory is 
in the symmetric phase with the minimum of the potential at p = 0. As the critical 
temperature is approached from above the mass term is given by 

m R {T) = U t (0,T) = ^^, (3.17) 

a result first obtained in [|J. In the language of critical exponents, the behaviour of m R 
corresponds to an exponent 

v = lim ;ri L — K J[, = 1. 3.18 
t^t c ,. d [ln(T - T cr )} y ' 

This result is in agreement with calculations of v in the context of the three-dimensional 
zero temperature theory Ul5 |. The quartic coupling in the symmetric phase is given by 



1 16vr 2 T 2 - T 2 

\ R (T) = U"(0,T) = -— 3 (3.19) 

It goes to zero as the critical temperature is approached with an exponent 

C = lim , ri L /rTi Hy J J N1 = 1. 3.20 
s T^T cr d [ln(T — T cr )] v ; 

The dimensionless ratio ^f^py approaches a temperature independent value at T cr 

llm Mill = (3.21) 
T^r cr m R (T) N ' K J 

a result which has an analogue in the language of three-dimensional field theory. (For 
an extensive discussion of this point we refer the reader to [10 , 15| |18| . ) It is interesting 



to recall at this point that the ratio ^ diverges near T cr and this is the reason for 
the breakdown of the standard perturbative calculations. We see how this divergence 
disappears when the bare quartic coupling is replaced by the renormalized one. At T cr 
the p 3 (0 6 ) term is the first non-zero one, with a coefficient given in ( |3.16|) . For T < T cr 
the potential has a minimum at a non-zero value 

Po(T) = ^(T 2 -T 2 ), (3.22) 
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and assumes the form 



U( P ,T) 



1 64tt 2 1 



(p-Po(T)f 



(3.23) 



after the addition of a p-independent constant. Notice that the minimum of the potential 
(the order parameter for the phase transition) approaches continuously zero as the critical 
temperature is approached from below (justifying our conclusion for a phase transition of 
the second kind). Expression ( |3.22|) predicts a critical exponent 



d 



Wp (T) 
lim — r^- , , 

t— >T cr d [ln(T cr - T)} 



0.5 



(3.24) 



again in agreement with the large N results of three-dimensional field theory. Another in- 
teresting point is the zero value of the second p-derivative of the potential at the minimum 
for T < T cr . This reflects the fact that the 4 theory is infrared free in the spontaneously 
broken phase due to the Goldstone modes. At non-zero temperature the dependence of the 
running coupling on the infrared cutoff k is linear [|l(J (compared to the logarithmic run- 
ning at zero temperature). This results in a zero renormalized coupling Xr{T) = U"(po, T) 
in the limit k — > 0. We see how this physical effect is naturally incorporated in our result 

A remark is due at this point about the form of the potential ( |3.23| ), since, obviously, 
the requirement of convexity is not satisfied. The paradox is resolved if one realizes that 
our discussion, and therefore the result (|3.23|), is valid only in the "outer" region of the 



potential p > p . For p < p the true saddle point of the functional integral ( |2.10[ ), ( |2.11| ) 
corresponds to C = and B given by the minimum of V(B). This saddle point is closely 
related to the spin wave solution of ref. [ HI , [12]. It results in the convex envelope of the 
non-convex part of (|3.23| ). 

The calculation of terms with higher powers of p in the effective potential requires 
additional terms in the high temperature expansion. Inserting (|3.6| ) and (|3.?1 ) into (|3.5|) 
and using ( |3.11|) , (|3.12| ) for the definition of p 2 R and \r, one obtains 



U'( P ,T) 



1 N 
32^2 



In 



M 2 



(dT) 2 + U>(p,T), 







T 2 
N — 


Pr 






24 


Xr_ 



+ p-N^y/U'(p,T). 



Differentiation with respect to p gives the p— and T-dependent quartic coupling 

Xr 



U"(p,T) 



1 



32tt 2 



ln(- { 



M 2 



U'(P,T) 



+ 



2nT 



(cT)*+U'(p,T)) (cT)Z+U>(p,T) ^ ^W(^F) 



(3.25) 
^3.26) 



It is interesting to study its dependence on the effective infrared cutoff U'(p,T). With 



| In U' one finds 







dt 



-M' 1 



N 
16^2 



n\2 



U\U' + 2(5T) 2 
JjJ> + (cT) 2 ) 2 



+ 7T- 



T 



(3.27) 
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We observe the transition from the usual logarithmic four-dimensional running | for U' ^> 
(cT) 2 to the effective linear three-dimensional running |H| |19[ for U' <C (cT) 2 . 



The critical temperature is again given by ( |3.14| ). For T near T cr and small p we can 
solve ( |3.25| ) approximately. Including the next to leading term we find, with M = cT, 

(3.28) 

with higher order terms involving larger powers of n ^t 2 ( T ~^ cr + ff)- The second term 
in ( 3.28 ) introduces corrections to the the effective potential ( |3.16|) which are not universal 
with respect to the three-dimensional effective theory. These corrections are sensitive to 
the properties of the full four- dimensional theory: they depend on \r. We note that the 
range where the pure three-dimensional behaviour dominates 

, n n 24 , W\ R T 2 

|T 2 - T 2 + — p| < o * (3.29) 



N 



is proportional to \r and, thus, very small for small \r. In the symmetric phase (T > T cr ) 
and in the critical region T — > T cr , the corrections to the first three terms in the p 
expansion are suppressed by additional powers of — jtr^- Keeping the leading contibution 
we obtain for the p 4 (0 8 ) coupling 

^1(0,7) = -^. (3.30) 

The negative sign of the s term does not imply the presence of a second minimum. It is 
obvious from eq. (|3.25 ) that above T cr there is only one minimum at p = 0. For T < T cr 
the minimum of the effective potential is given by ( |3.22| ) and the potential assumes the 
form 

V(P, T) = j~ 2 ^f^ (P - Po(T)) 3 - ^^20487^ (p - p (T)) 4 . (3.31) 



Eqs. ([3T28D , (ggg ) can be directly compared to ( gT6|) , (|3^ ). We conclude that the 



universal behaviour near the critical temperature (which determines the nature of the 
transition and the critical behaviour of the system) is effectively three-dimensional and 
is given by the first three derivatives of the potential around its minimum (the "rele- 
vant" terms). Higher derivatives of the potential do not affect the universal behaviour. 
They simply involve secondary non-universal corrections to the effective potential. One 
should keep in mind, however, that the non-universal corrections (which reflect the four- 
dimensional nature of the theory) must be taken into account for temperatures away from 
the critical one. 

In summary, the universal part of the effective potential for the 4 theory in the critical 
region is given by eq. ( 3.16Q for large N. In the following section, we shall reproduce this 
result by solving the Schwinger-Dyson equations for the high temperature theory near the 
phase transition. 



§Note that, for U' 3> (cT) 2 , \fJt%\, one has approximately U' ~ Xrp and ( 3.27 ) reproduces the result 
of in leading order. 
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4. The Schwinger-Dyson equations 



In this section we are interested in verifying the results of our previous discussion through 
a method more closely connected to the standard perturbative evaluation of the effective 
potential J7|, || ^ [|]. It is well known J3J that the Schwinger-Dyson equation for the 
propagator, for a 4 theory in the large N limit, can be reproduced by the summation of an 
infinite class of diagrams, usually characterized as "daisy" or "cactus" graphs [f|, |]. The 
solution of this equation gives the value of the critical temperature and the temperature 
dependence of the mass term in the critical region. We have seen in the previous section 
that the first three terms of the effective potential in the critical region (|3.16| ) do not 
involve the parameters of the four- dimensional theory, but are fully determined by the 
dynamics of the effective three-dimensional theory. Moreover, the form of the potential 
predicts a characteristic temperature dependence for its second and third p-derivatives 
(which correspond to the renormalized four and six point vertices). These results have so 
far no analogue in the context of perturbation theory at non-zero temperature, since the 
renormalization of the quartic coupling has not yet been considered [|], 3 . It is interesting 
to investigate whether the behaviour predicted by eq. (|3.16|) can be reproduced through 
the study of the relevant Schwinger-Dyson equations, since a solution of these equations 
automatically accounts for the renormalization of the couplings. 

In fig. 1 we give a graphical representation of the first three Schwinger-Dyson equations 
for the O(N) <ft A theory in the symmetric phase for large N. (For a derivation of the 
Schwinger-Dyson equations in the functional integral formalism see ||20|| .) Only the leading 
terms in 1/N have been included. (We remind the reader that the vertices come with an 
power of N given by ( |2.13| ).) The left-hand-side of the equations involves the first three 
p-derivatives of the potential at zero external momentum (the renormalized mass and four 
and six point vertices). Simple vertices denote bare quantities, and lines with a circle full 
propagators. The representation is schematic and the signs and combinatoric factors for 
the various terms have not been included. The incorporation of non-zero temperature 
effects in the formalism proceeds exactly as in the previous section. By recalling eqs. 
(|3.3|) , ( |3.4[) we can express eq. (a) of fig. 1 as 



U'(0,T) 



y A/(T) - f 



+ ^X[I 1 (U'(0 : T),T)-I(T)]. (4.1) 



For the solution of the above equation we proceed as for the solution of (|3.5| ): We use 
eqs. (|3.6| ) and ( |3.7|) , introduce the renormalized quantities \r according to (|3.11| ), 
( |3.12j ), define the critical temperature (|3.14| ), and finally obtain the standard result [[| 

"ttT 2 -T 2 



U'(0,T)=m 2 R (T) 
Eq. (b) of fig. 1 reads in our notation 



- cr 



3 T 



4.2 



N- 

U"(0,T) = A - -XU"(0,T) I 2 (U'(0,T),T), (4.3) 
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where I2 was defined according to ( |3.4j ). Making use of the high temperature expansion 
(|OD we find 

U"(0, T) = — - — r . (4.4) 

1 + £\T[U'(0,T)}^ 

In the critical region T — > T cr the second term in the denominator of the last expression 
(which comes entirely from the effective three-dimensional theory) diverges. Substitution 
of ( ^4.2| ) gives 

1 16vr 2 T 2 - T 2 

U"(0,T) = X R (T) = - — (4.5) 

Notice that all reference to A has again disappeared in the above equation. In a similar 
way we obtain for eq. (c) of fig. 1 

N - 1 
U"'(0, T) = -—\U'"(0, T) I 2 (U'(0, T), T) + NX [U"(0, T)] I 3 (U'(0, T), T). (4.6) 

Through use of Q, (Q, Q, Q we find 

U"\0,T) = ±- 2 l28n 2 ± (4.7) 

Reconstruction of the effective potential in the symmetric phase, through a Taylor 
expansion around the origin, results in eq. ( |3.16|) . It should also be noted that eqs. 



|4.1|) , ( f4.3|) , (|4.6|) can be obtained from fl3.2p by differentiation with respect to p. We 



have, therefore, shown that the solution of the Schwinger-Dyson equations in the large N 
limit verifies the results of the discussion of the previous section. Moreover, it provides 
helpful insight on how the first three terms of the effective potential for T — > T cr become 
independent of the four- dimensional parameters, due to the singular infrared behaviour of 
the effective three-dimensional theory. In the following section we shall use the evolution 
equation for the effective potential in order to derive results beyond the leading order in 
1/N. 



5. Beyond leading order in 1/N 

The picture for the high temperature phase transition for the 0(iV)-symmetric scalar 
theory is completely satisfactory for large values of N. In order to deal with realistic 
theories, e.g. N = 4 for the scalar sector of the electroweak standard model or the QCD 
phase transition for vanishing up and down quark mass 0, one has to go beyond the 
leading order approximation and include corrections ~ 1/N. Near the phase transition 
the three-dimensional behaviour is dominant and any useful method should account for a 
proper treatment of the three-dimensional 0(iV)-symmetric scalar theory. In particular, 
one expects the critical exponents is, ft, ( to deviate from their simple values ( |3.18|) , (|3.24| ), 
( |3.2C| ), so that the effective potential becomes non-analytic in T — T cr . 

Going beyond the leading result in 1/N is not as straightforward as it may seem at first 
sight. The difficulty lies in the infrared behaviour of the three-dimensional theory. At the 
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critical temperature, this leads to the infrared divergencies of standard high temperature 
perturbation theory. A systematic expansion of (|2.10 ), ( 2.11|) beyond the saddle point 



approximation is not easy, even for k = 0,P(— □) = — □ and (j) a = const. As far as the 
universal behaviour at the critical temperature is concerned, one can rely, of course, on 



results of the purely three-dimensional 1/N expansion [TIJ. This should correctly predict 
the first 1/N correction to the critical exponents. However, for an estimate of either the 
critical temperature, or the behaviour of U(p, T) away from the critical temperature, or 
the proportionality constants in relations such as tur ~ (T — T cr ) v ', the properties of the 
underlying four-dimensional theory play a role. For these questions the four-dimensional 
1/N expansion at no n- vanishing temperature is needed. The same remark applies to other 
standard methods employed in three-dimensional field theory such as the 4 — e or 2 + e 
expansions or the strong coupling results on the lattice. They give reliable results for the 
critical behaviour but cannot be employed for the behaviour away from T cr , where the 
matching with the four-dimensional behaviour is needed. 

One could, in principle, consider the Schwinger- Dyson equations and include terms 
subleading in 1/N. These terms, however, turn the Schwinger-Dyson equations into 
a complicated coupled system, as illustrated in fig. 2. The computation of the two 
point function involves now the full four point vertex and this difficulty propagates to all 
vertices. As we have seen before, a replacement of the full four point vertex by the bare 
coupling A does not reflect properly the infrared behaviour near T cr , where the appropriate 
renormalized coupling Xr(T) vanishes. More precisely, the last diagram in fig. 2 with A 
instead of the full four point vertex is infrared divergent at T = T cr , whereas the correct 
expression is infrared finite. 



All these effects are naturally incorporated in the evolution equation ( 2.27Q , which 
has a simple generalization O, [TI| for arbitrary N. When wavefunction renormalization 
is neglected, the generalized evolution equation reads 

d TT , 1 f d d q ( N-l 1 \,„ d nf „.2x 



k dk Uk{p) 2 h (2tt)<* [Pig*) + U' k (p) + P{f) + U>(p) + 2U»(p)p) k dk P{q 

f A2 a --i I N ~ 1 1 ^ u d 

- Vd Jo ^ [pjx) + U' k {p) + P{x) + U' k {p) + 2U» k (p)p k dk P{x){5A) 



where Vd is given by ( |2.2£ ). One can easily identify in the above expression the con- 



tributions of the N — 1 Goldstone modes and the radial mode. In the limit iV — » oo 
the contribution of the Goldstone modes dominates and the evolution equation becomes 



identical to Q2.27|) . The generalization of the last equation in order to include non-zero 
temperature effects is straightforward and follows the discussion in the end of section 2. 
For the four- dimensional theory one finds 



^ This equation can also be interpreted as a truncation of an exact evolution equation pl[ (for arbitrary 
N) which is related to an improved average action p2| . 
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N - 1 

+ 



\P(cp + 4vr 2 m 2 T 2 ) + U' k (p, T) P{q> + 4vr 2 m 2 T 2 ) + U' k (p, T) + 2U k \p, T)p J 

(5.2) 

For k ^> T the summation over the discrete values of m in the above equation is equal to 
an integration over a continuous range of q°, up to exponentially small corrections. In this 
limit eq. ( |5.2| ) is identical to ( |5.1|) with d = 4. In the opposite limit k -C T the summation 
over m is dominated by the m = contribution and the evolution becomes effectively 
three-dimensional. The full solution of (|5.2|) interpolates between the two regimes and, 



thus, matches the four- dimensional behaviour with the three-dimensional one. 

For the study of the phase transition we are interested in solving ( |5.2|) near the min- 
imum of the potential for different temperatures T. In order to do so we parametrize 
Uk(p,T) in terms of the location of its minimum p (k,T) and its successive p-derivatives 
at the minimum U k (p ,T),U k (po,T)... This turns the partial differential equation (|5.2|) 
into an infinite system of coupled differential equations, which can be solved approximately 
by truncation at a finite number of equations. We have seen in the previous sections that 
the order of the phase transition as well as the critical behaviour of the system are deter- 
mined by the "relevant" first three p-derivatives of the effective potential. It is sufficient, 
therefore, to keep the first three coupled differential equations of the infinite system. We 
give these equations without presenting the details of their derivation here, since they 
have been repeatedly derived elsewhere |TT|, [II], . In the spontaneously broken regime 
(po 7^ 0, U' k — 0) the equations for p (k,T), U' k r (po,T) and U k (po,T) read 

^ = - (3 + ^) LfrUfa, T) + (N- l)Li(0, T) } (5.3) 

dJ § = ~ v*\Uff{ (3 + ^) 2 Lt(2U' k W T) + (N- 1)^(0, T) } 

+ v 4 k* (V - ^f^) iAWito T) (5.4) 
d -^f = 2„ 4 ATW j (3 + ^) 3 ^(2^p ,T) + (N — 1)^(0, T)} 

- to<UZU?L + L&U'^T) + (N- 1)L*(0, T) J. (5.5) 

In the symmetric regime (po = 0) the corresponding equations for U' k (0,T), U k (0,T) and 
Ujl'(0, T) read 

dlV 

^= v 4 k\N + 2)U>>Lt(U>,T) (5.6) 



16 



dm 



dt 
dUL 



1 - - v,(N + 8) [UfflAW* T) + v 4 k\N + 4)U£Li(U' k , T) (5.7) 
2v,k- 2 (N + 2Q)[Un 3 Li(UiT)-3v i (N+U)U'^'Li(UiT). (5.8) 



dt 

The dimensionless integrals L^(w,T) 

Li(w, T) = -nk 2n -^T £ J d 3 q(P(<f + 4vr 2 m 2 T 2 ) + w) " (n+1) k-^P(<f + 4vr 2 m 2 T 2 ) 

(5.9) 

have been extensively discussed in |]10|| . (Notice also that these integrals are proportional 
to the logarithmic derivative with respect to k of the integrals I n discussed in section 3, 
with q 2 replaced by P{q 2 ) in the propagator.) They exhibit a behaviour as a function 
of j which interpolates between being purely four-dimensional for T <^ k and effectively 
three-dimensional for T ^> k. This induces a continuous transition from four to three for 
the effective dimensionality of the evolution equations. 

The high temperature phase transition for the 0(iV)-symmetric <f) A theory was inves- 
tigated in detail in ref. |i~0| . We shall not repeat the technical discussion here and we 



refer the reader to [JTOj for a thorough presentation. We briefly summarize the main re- 
sults in order to point out the consistency with the conclusions of the previous sections. 
The evolution equation was studied through the parametrization of the average poten- 
tial in terms of its minimum and its derivatives at the minimum, as explained earlier in 
this section. The resulting infinite system of differential equations was truncated after 
two equations and only the first two derivatives of the potential were considered. This 
amounts to solving the equations ( |5.3| ), (|5.4j ) and ( |5.6| ), ( |5.7| ) without the terms involving 



U'k'. The evolution equations were solved at various temperatures in order to obtain the 
parameters of the theory in the limit k — > 0. In this limit the average potential becomes 
equal to the effective potential, and its minimum gives the ground state of the theory and 
its derivatives the renormalized mass and couplings. The phase transition was studied 
by defining the renormalized theory in the spontaneously broken phase at T = 0, and 
following the evolution of the minimum of the effective potential for increasing tempera- 
ture. The phase transition was found to be of the second order with a critical temperature 
given by (|3.14j ) for small values of the quartic coupling, independent of N. The theory 
was shown to have a vanishing quartic coupling at and below the critical temperature for 
any N ^ 1, in agreement with ( |3.23|) . The universal critical behaviour was found to be 
determined by a fixed point corresponding to a second order phase transition. The same 
fixed point determines the critical behaviour of the effective three-dimensional theory, in 
exact agreement with the discussion in sections 3 and 4. An extrapolation of the numer- 



ical results for the critical exponents in [10|, to large values of N, indicates agreement 



with ( |3.18| ), ( |3.20| ), (|3.24j ). Evaluation of the analytic expressions for the exponents in 
[|T0f for iV — > oo reproduces the same values. Finally the numerical values for the ratio 
\^{Tcr)Tcr consistent with (|3.21|) for large values of N. 



We have repeated the calculation of the critical exponents by considering the first three 
derivatives of the potential. This was done by solving the system of evolution equations 



17 



( |5.3| )-( |5T5D and ( |5.6| )-( |5T8| ). In this way all the parameters "relevant" for the universal 
critical behaviour (apart from the wavefunction renormalization which is non-zero for 
finite N) were taken into account. The results of this calculation for the exponent (3 (for a 
value 6 = 3, see (|2.4Q ) are presented in table 1. The results for the other exponents satisfy 
v = 2(3 (in agreement with the scaling law for zero wavefunction renormalization) and ( = 
v. In the same table we list values for the exponent (3 obtained through the 1/N expansion 
to order 1/N 2 for the three-dimensional theory |23|| . We observe exact agreement for large 
values of N, which becomes bad for small values of N. For N = 3 we also give the value 
for the same critical exponent obtained through other methods (e-expansion, summed 
perturbation theory in three dimensions, lattice calculations) Hl5| . It is clear that our 
result is in good agreement with this value, while the 1/N expansion fails for small N. The 
small remaining discrepancy is due to the omission of the wavefunction renormalization 
from the evolution equations. The study of the wavefunction renormalization effects is 



under way J24 



We summarize this section by concluding that the solution of the evolution equation 
fully confirms the results of the previous two sections, while it provides the means for the 
detailed study of the phase transition for small values of N. We also emphasize that this 
method gives a reliable description for temperatures away from the critical one [ELO] . 



6. Conclusions 

We have studied the high temperature phase transition for the iV-component </> 4 theory 



in the large N limit, in order to confirm the results of ref. |TIJ in which the method of 
the average potential was employed. The investigation was carried out through use of 
different methods, such as the standard large N evaluation of the functional integral for 
the effective potential and the solution of the Schwinger-Dyson equations. Our results 
fully confirm the results of [|H]]. The phase transition is of the second kind and falls in the 
same universality class as the phase transition of the three-dimensional 0(iV)-symmetric 
a-models. The critical behaviour is determined by the first three derivatives with respect 
to 4> 2 of the effective potential around its minimum, and it is effectively three dimensional. 
Explicit expressions have been obtained for the critical temperature and the effective 
potential in the critical region. As the critical temperature is approached from above, the 
renormalized mass term m^(T) in the potential approaches zero with a critical exponent 
v = 1. The renormalized quartic coupling \r(T) approaches zero also, with a critical 
exponent £ = 1. The ratio ^^py remains finite and, as a result, no infrared divergencies 
arise. At the critical temperature the first non-zero term in the potential is the 6 term, 
whose coefficient has been calculated. Below the critical temperature the minimum of the 
potential (which gives the ground state and, therefore, the order parameter of the theory) 
moves continuously away from zero, with a temperature dependence characterized by an 
exponent (3 = 0.5. All the above is fully consistent with the results of ref. Jl0| and 
three-dimensional field theory. 

Going beyond the leading 1/N result, in order to discuss small, physically relevant 
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values of N, is more difficult and the renormalization group becomes an indispensable tool. 
We have demonstrated the power of the method of the average potential by calculating 
the critical exponents for various values of N. We have found perfect agreement with 
calculations using the 1/N expansion in three dimensions. Morevoer, we have shown that 
the average potential gives the correct answer even for small values of N, where the 1/N 
expansion fails. The final ingredient for the complete description of the phase transition 
is the wave function renormalization, and this investigation is under way 
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Tables 



N 




1 nn 


U.OU 


n ^nW 

U.OLT ' 


20 


0.48 


0.48^ 


10 


0.46 


0.46^ 


4 


0.41 


0.34W 


3 


0.40 


0.25^ 0.38^ 



Table 1: Critical exponent (5 for various N (b = 3). For comparison we have listed the 
results of calculations of the same exponent in three-dimensional field theory: 

a) From 1/N expansion to order 1/iV 2 fl23 |. 

b) Average value of results from e-expansion, summed perturbation theory in three di- 



mensions and lattice calculations 15 
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Figures 



Fig. 1 Schwinger- Dyson equations in the large N limit. 



Fig. 2 Schwinger- Dyson equation for the two point function. 



23 



